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The vibrational modes of Jahn- Teller molecules are affected by a Berry phase that is associated 
with a conical intersection of the adiabatic potentials. We investigate theoretically the effect of this 
Berry phase on the electronic transport properties of a single E®e Jahn- Teller molecule when the 
tunneling electrons continually switch the molecule between a symmetric and a Jahn- Teller distorted 
charge state. We find that the Berry phase, in conjunction with a spectral trapping mechanism, 
leads to a current-blockade even in regions outside the Coulomb blockade. The blockade is strongly 
asymmetric in the gate voltage and induces pronounced negative differential conductance. 
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I. INTRODUCTION 

The vision of molecular electronics has stimulated 
great interest in both the experimentaUii^i^ii^i^ and the 
theoretica&^i^iiSiiiiii^ii^ii^ understanding of electronic 
transport through single molecule devices. The cou- 
pling of electronic degrees of freedom to few, well- 
defined molecular vibrations, a property that discrim- 
inates transport through single molecules from that 
through other nanostructures such as quantum dots, sets 
the stage for observing novel quantum-transport phe- 
nomena. ^5, Qjjg Qf ^i^g simplest and most intensively stud- 
ied models of such a device is a molecule with a single 
electronic level and a one-dimensional potential energy 
surface for the nuclear displacements. This model al- 
ready gives rise to a large variety of interesting trans- 
port phenomena such as vibrational sidebandsjii^ elec- 
tron shuttles,^ Franck-Condon blockade,— avalanche-like 
transport, pair tunneling, and dynamical symmetry 
breakingfi^ Generalizations of this model to two or more 
electronic levels have been suggested in order to account 
for non-degenerate but competing molecular statesi^ and 
for degenerate electronic states on a molecular dimerJ^ 

Extensions of the simple molecular model to higher- 
dimensional potential energy surfaces may seem to 
be straightforward without yielding qualitatively new 
physics. This is, however, not the case for molecules with 
symmetry-induced degeneracies of both electronic and 
vibrational states. For such molecules, the Jahn- Teller 
effecti^- leads to nontrivial vibrational dynamics, origi- 
nating in the conical intersections of adiabatic potentials 
and their associated geometric (Berry) phases. We ex- 
pect the transport properties to be particularly sensitive 
to this intriguing dynamics when the tunneling of elec- 
trons between molecule and leads switches the molecule 
between a symmetric and a Jahn- Teller distorted charge 
state. 

In this paper, we investigate the consequences of such 
a Berry phase for transport through single Jahn- Teller 
molecules. We focus on the E®e Jahn- Teller effect, for 
it is the simplest example with a conical intersection of 
the potential energy surfaces. This type of Jahn- Teller 
effect occurs for a large variety of triangular (X3), tetra- 
hedral (ML4 — a transition metal ion M surrounded by 




FIG. 1: Top: Adiabatic potentials of the E®e Jahn- Teller 
distorted (charged) molecule (left) and the undistorted (neu- 
tral) molecule (right) as a function of {Qg,Qe). Bottom: 
Displacement patterns of the degenerate 9 and e vibrations 
for octahedral MLg molecules. 



a tetrahedron of ligands L), and octahedral molecules 
or complexes MLg. We show in section III Bl that the 
Berry phase induces a nontrivial selection rule for transi- 
tions between vibronic excitations of different molecular 
charge states. This selection rule, in combination with 
generic properties of the energy spectrum of the Jahn- 
Teller molecule, leads to the formation of trapping states 
in the sequential tunneling limit (Section IIIip . The con- 
sequences are a current-blockade outside the Coulomb 
blockade window, asymmetry of the current-voltage char- 
acteristics, and strong negative differential conductance 
(Section IIV[) . The appearance of trapping states also 
strengthens the role of higher order processes, which in- 
troduce a second time scale that is much larger than the 
scale on which sequential tunneling takes place. The in- 
terplay of these two scales renders the details of the tun- 
neling dynamics irrelevant and has a non-negligible effect 
on the line shape of the differential conductance peaks. 
We conclude the paper with a discussion of some exper- 
imental issues (Section |V| and an outlook to interesting 
problems related to electronic transport through Jahn- 
Teller molecules (Section IVI[) . 
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II. JAHN-TELLER MOLECULE 
A. Model 

The E®e Jahn- Teller effect generally occurs in 
molecules where two degenerate electronic states are cou- 
pled to two degenerate vibrational modes 9 and e. In the 
case of an octahedral MLg molecule, the relevant elec- 
tronic states are of dx^-y'^- and d^z'^-r^-ty^e, and the 
degenerate vibrations correspond to the symmetric oscil- 
lations of the ligands as they are pictured in Fig [1] For 
definiteness, we assume that without an excess charge 
on the molecule the degenerate electronic E states are 
empty 1^ With this choice, the molecular Hamiltonian 
isii 

HE<^e =eV^n^ + ^ {P^ + Ql + P'; + Ql) 

+ Xhhj [Qedx + QeCTy) rid- 

We use a pseudospin notation with Pauli matrices <7i for 
the electronic states I +), |— ). The dimensionless strength 
of the linear coupling of electronic and vibrational de- 
grees of freedom is A, and lo is the frequency of the vibra- 
tional mode with effective mass M and oscillator length 
^osc = {h/AIwY^'^ . The dimensionless coordinates Q and 
momenta P of the normal modes 9 and e are measured 
in units of ^osc- The gate voltage Vg tunes the energy 
of the degenerate electronic levels. = means that 
for vanishing bias, the molecular levels are degenerate 
with the Fermi levels of the electrodes. Electronic oc- 
cupations different from = (neutral molecule) and 
— I (charged molecule) are effectively projected out 
by assuming a large charging energy U . 

The potential surfaces of both charge states are U{1) 
symmetric in the {Qe,Qe) plane and thus have a con- 
served angular momentum quantum number. We shall 
refer to this angular momentum as a pseudo angular mo- 
mentum to make clear that it does not describe rota- 
tions in real space but rather a continuous change of the 
shape of the molecule's distortion. A vibronic state of 
the neutral molecule is labeled by its integer pseudo an- 
gular momentum I and radial excitation nr. As it is well- 
known from the occupation number representation of the 
two-dimensional harmonic oscillator, its energy depends 
on the total number of vibrational quanta no matter in 
which direction of the oscillator they are excited. The 
energy manifold E — {N + l)huj is iV-fold degenerate. In 
the pseudo angular momentum representation, the de- 
generacy is reflected in the relation = Ei^2,n^-i, 
compare Fig. [5] (a). 

In the charged configuration, on the contrary, with an 
explicit coupling of vibrations to the pseudospin, it is the 
total pseudo angular momentum j — I + ^cr that is con- 
served. That the conserved pseudo angular momentum 
must be a half-integer quantity can also be seen by in- 
spection of the adiabatic potentials, which are computed 
by diagonalizing the potential energy of He^c for fixed 
nuclear displacements {Qe,Qe)- The two sheets of the 
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FIG. 2: Vibrational spectra (a) of the neutral and (b) of 
the charged molecule vs. Pseudo angular momentum. Only 
energies for positive pseudo angular momenta are shown as 
the spectra are symmetric under the reflection / t-^ —I. 

— 1 adiabatic potential have, in contrast to the single 
sheet of the neutral state, a conical intersection at the 
origin. Along any loop enclosing the conical intersection, 
the electronic pseudospin of the adiabatic eigenstates ro- 
tates by 27rji^ The adiabatic eigenstates consequently ac- 
quire a Berry phase of tt along such closed loops, which 
must be offset by the vibrational wavefunction in order 
to make the full, electronic plus vibrational, wavefunc- 
tion single valued. This implies antiperiodic boundary 
conditions for the vibrational wavefunction and thus half- 
integer pseudo angular momentum j. The spectrum of 
the charged configuration is plotted in Fig.[2](b) for an in- 
termediate value of the electron-phonon coupling, A = 1. 
In contrast to the regularity of the oscillator spectrum, 
the Jahn- Teller effect lifts the large degeneracy; each 
state is only doubly degenerate, Ej^n^ = E-j.n, ■ The 
Jahn- Teller effect also lowers the ground state energies 
of each pseudo angular momentum column with respect 
to the harmonic oscillator energies, which is one of the 
keys to the current-blockade mechanism discussed in this 
paper. 

The molecular system is coupled to two non-interacting 
Fermi-gas electrodes labeled by a = L, R. We define elec- 
tron annihilation operators Caks for a particle in lead a 
with wavevector k and spin s, and d±^s for a particle 
on the molecule. The tunneling between the molecule 
and the electrodes is described by a conventional tun- 
nel Hamiltonian, Ht = taJ2a.k,s'^aks{d\_ ,, -\- J -I- h.c. 
Jahn- Teller distortions usually occur on a sub- Angstrom 
scale^i which is small compared with the scale on which 
electronic wavefunctions fall off. For this reason, we can 
neglect the dependence of to on the vibrational coordi- 
nates. The full Hamiltonian of the single molecule device 
is 

H = He^c + ffloads + ^^T, (2) 

where iJioads = Y,ks ^kclks^aks IS the Hamiltonian 
of the electrodes. The states of the molecule will be la- 
beled |n.d; j, nr), rid is the charge on the molecule, j is the 
pseudo angular momentum {I in the neutral state), and 
fir enumerates the radial excitations. Our discussion is 
restricted to the weak tunneling case, where the tunnel- 
ing induced broadening of the molecular states F is small 
compared with the temperature T, and the predominant 
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transport process is sequential tunneling of electrons. 
We will only consider a regime of intermediate electron 
phonon coupling, where the system is not in the Franck- 
Condon blockade. The Franck-Condon blockade effect 
would not alter the validity of our conclusions but merely 
complicate the discussion unnecessarily. 



B. Selection Rule 

Any computation of transport properties via an ex- 
pansion of the system's reduced density matrix in pow- 
ers of the tunnel Hamiltonian requires the computation 
of the matrix elements of Ht- An expansion of the 
Jahn- Teller states |l;j, rir) in terms of harmonic oscil- 
lator eigenstates \l,n'^) times the two-dimensional elec- 
tronic manifold |±) shows that the electron-phonon inter- 
action due to the E(^e Jahn- Teller effect only couples the 
states with the states + of the non- 

interacting systemji^ The associated half-integer pseudo 
angular momentum number of the Jahn- Teller state is 
j — / + 5. The eigenstates of the Jahn- Teller active 
charge state thus have the general form 



|1; j,"-r) = 



^ri+)b-^'"r)+i?ri-)ij+j.^r) 



(3) 

for suitable expansion coefhcients A-"' and B~^\ The 
matrix elements (0; Z, nj. [ \ 1; j, n^) will therefore van- 
ish unless the pseudo angular momenta fulfil 
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(4) 



This condition causes many zeroes in the tunnel matrix 
and has a great influence on the transport properties of 
the device as will be shown in the following sections. 
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FIG. 3: Differential conductance vs. gate voltage Vg and 
source-drain voltage Vsd for A = 1, ksT = 0.01 hu. (a) 
T = 10 TO, (b) T — 5000 To. White correspond to positive, 
Black to negative values. Neutral Gray is zero differential 
conductance. 



The implications of this effect for the differential con- 
ductance are the formation of structures like those in 
Fig. [21 a suppression of the stationary current outside 
the Coulomb blockaded region. The suppression, how- 
ever, is not complete; a small stationary current is admit- 
ted due to the presence of higher order processes such as 
co-tunneling, relaxation to thermal equilibrium, or mix- 
ing of states by external perturbations, which effectively 
break the strict selection rule, allow to release the sys- 
tem to a vibrational state that can then be left by tun- 
neling. Since these processes are of higher order in to, 
they define a second time scale r. This time scale, which 
describes the time that the system will effectively spend 
in a trapping state, is much larger than the time scale of 
sequential tunneling, given by Ht, Tq^ := 27ri^iQ, where 
1/ is the density of states in the electrodes. 



III. CURRENT BLOCKADE 

The main result of our paper is a current-blockade ef- 
fect outside the Coulomb blockade diamonds. In certain 
areas of the (Vg, Vsd) plane, the system exhibits "trapping 
states" that, if populated, lock an excess charge on the 
molecule and, due to the large charging energy, obstruct 
any stationary current to flow through the device. These 
states are vibrational eigenstates of the Jahn- Teller active 
molecule with pseudo angular momentum \j \ > i. The 
effect becomes more pronounced for larger pseudo angu- 
lar momenta. For certain voltage parameters that lie well 
outside the Coulomb-blockaded region, an electron tun- 
neling from the source electrode can excite the molecule 
into such states, but it cannot tunnel out towards the 
drain electrode. The trapping effect itself comes about, 
because these states become isolated in the tunneling dy- 
namics: either a transition to a neutral target state is 
forbidden by the selection rule or it is forbidden by con- 
servation of energy. 



A. Formalism 

The technical and formal discussion of the trapping ef- 
fect is subtle, because due to the high level of degeneracy, 
one has to employ full Master equation dynamics instead 
of intuitive rate equations for sequential tunneling. The 
qualitative picture can nonetheless be grasped by rate 
equations in the pseudo angular momentum basis. 

We confine ourselves to the sequential tunneling 
regime, which means that transitions will only take 
place between different charge states. We use an expan- 
sion of the reduced density matrix p in the tunneling 
Hamiltonian, * employ both the Markov and the Born 
approximation, and eventually arrive at a Liouville-von 
Neumann equation, p = — «[i?ioads + He^c, p] + C,p- The 
operator C has the well-known Lindblad form for dis- 
sipative linear operators and describes the tunneling in 
second order perturbation theoryiJ^ Since we are only 
interested in the stationary transport properties, we can 
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safely neglect those off-diagonal matrix elements of p that 
belong to states of different energies and thus effectively 
produce a rate equation dynamics between different en- 
ergy manifolds. The stationary von Neuman equation 
reduces with this approximation to Cp = 0. But since 
the system under consideration has degenerate states, 
the equations can, at this point, not be further simpli- 
fied to rate equations for the states |l;j,ni) and |0;Z,nJ.). 
We do, however, know that the density matrix is a self- 
adjoint operator and hence admits an eigenbasis in the 
degenerate subspaces for each set of voltage parameters 
(V^d,Vg). In this basis, the von Neumann equation as- 
sumes the form of a rate equation, where the rate for a 
transition between two eigenstates \P) and I7) is given by 
a golden-rule expression, 

MP-'' ■.^'^vtl\{^\dl + dl\[3)\\ (5) 

which is of order 0{H^), multiplied by a Fermi func- 
tion that accounts for the availability of an electron or a 
hole in the leads. If we want to compute the rates be- 
tween different energy manifolds in the tunneling prob- 
lem, we shall consider the, unknown, diagonal basis of 
the degenerate subspaces of p at fixed voltages. Let 
7) — '^dn^^i \l,n-r) be one of these states in the neu- 
tral configuration and — a.j \ j,n[) + a_j |— j, rtj.) one 
in the charged configuration. The Franck-Condon factor 
is 

|(7|i?T|/3)P = | 5] <^,AK,Z|iJT|<,j)P 

(6) 

The selection rule, which has been shown to be strict 
in the pscudo angular momentum basis, has thus impli- 
cations on the transition matrix elements whatever the 
eigenbasis of p is. If the angular momentum of a charged 
state 1/3) is too large compared with those present in a 
neutral state [7), a transition may not be allowed. We 
will show in the next section that such forbidden tran- 
sitions between energy manifolds are the cause of the 
current-blockade that we observe analytically and numer- 
ically. 

B. Trapping states 

Considering a sufficiently large bias and a gate volt- 
age that places the ground state of the charged molecule 
slightly above the Fermi level of the drain electrode, the 
Fermi functions for relevant excitations are 1 — /source — 
/drain ~ 0, and wc Only need to discuss the rates for 
an electron to leave the molecule towards the drain elec- 
trode. In this configuration, the system is not in the 
Coulomb blockaded region and ought to be conducting. 
This is, however, not the case due to the selection rule 



Eq. ([4]) and the spectral structure of the molecule. Con- 
sider Fig. m in which the vibronic spectra of both charge 
states have been overlaid and shifted by eVg and ±iel4d 
along the energy axis to account for the electronic en- 
ergies of the charged and the neutral state respectively. 
Conservation of energy then implies that transitions are 
energetically possible if the charged target state lies be- 
low the neutral initial state in Fig. 2] (a) and vice versa in 
Fig. |4] (b), see the figure caption for more details. Take 
for example the state \(3) := |1;|,0), the charged state 
with pseudo angular momentum J = f and radial ground 
state. Once the system occupies this state, it will stay 
there forever; the electron cannot tunnel to the drain elec- 
trode, see Figs. Eland [6l Conservation of energy requires 
that any target state has to have less energy than so 
only / < 1 and suitable n are possible. The selection rule 
Eq. ([4]) only allows the dynamics to connect neighboring 
columns of Fig. [H But then all transition rates are zero. 
The only states that have non-vanishing Franck-Condon 
factors belong to energy manifolds that due to the spe- 
cific choice of voltages are out of reach energetically. The 
state \f3) is called absorbing or "trapping state", that is 
to say there is no escape rate but yet a finite rate for 
populating the state. The latter of the two properties 
of 1/3) is easily verified with Fig. |4] (a), which shows the 
shift of the two spectra for tunneling processes from the 
source electrode. In this configuration there are enough 
states from which \(3) can be populated, for instance the 
E = Ahuj manifold. 



The energy threshold above which a trapping state can 
be reached from the vibronic ground-state of the neutral 
system marks the onset of the current-blockade in the 
dl/dV diagram. It can thus cause negative differential 
conductance peaks, like those in Fig. [3l if the system 
were conducting below the threshold. Positive differen- 
tial conductance peaks, on the contrary, can be observed 
for the energy threshold above which the trapping state 
acquires a finite escape rate and is no longer trapping. 
The current-blockade effect, as it has been discussed in 
the previous paragraph also produces a significant asym- 
metry of the dl/dV diagram with respect to Ig. In Fig.[5l 
we plot the energies of the charged and uncharged states 
for voltage parameters that place the molecule closely 
below the Fermi level of the source electrode, that is for 
the same Vg as before but with opposite sign. In both 
processes, tunneling onto and off the molecule, there is 
no state of the charged system that is isolated due to 
the selection rule or the spectral distortion. The prop- 
erty of the Jahn- Teller effect to lower the energy of the 
distorted states is the actual reason for this asymmetry. 
If our Jahn- Teller molecule were such that it was Jahn- 
Teller active without an excess charge and symmetric with 
an additional electron, the situation would be inverted, 
current-blockade for positive, conducting at negative gate 
voltage. 



FIG. 4: Angular-momentum resolved spectra of the Jahn- 
Teller molecule (A = 1) (diamonds) and the neutral molecule 
(circles). Radial excitations are grouped for each pseudo an- 
gular momentum. Panel (a) shows the situation for tunnel- 
ing at the source electrode, panel (b) the situation at the 
drain electrode. The chosen parameters are eVad = 2hu!, 
eVg = —O.bhiu. Consider panel (a). An electron can only 
tunnel from the source electrode and induce the vibronic 
transition \0;q) ll;?') if the difference in the electronic 
energy, |eV^d — eVg, is larger than the difference in the vi- 
bronic energy, Eiqi — Eoq. More precisely, the Fermi function 
/(e(Vg — 5 Vsd) + Eiqi — Eoq) has to be positive. By shifting 
the neutral spectrum by ^eVsd and the charged spectrum by 
eVg, this relation translates into the intuitive rule that in the 
diagram the target state of a tunnel process has to lie below 
the initial state. Energetically allowed transitions are there- 
fore directed downwards in the figure. A horizontal transition 
corresponds to the borderline case of an electron sitting ex- 
actly at the Fermi level of the electrode. That conservation of 
energy has the form of an inequality is due to the filled con- 
tinuum of states in the Fermi-gas electrode. Panel (b) shows 
the corresponding energetic situation for sequential tunneling 
processes at the drain electrode. The presence of trapping 
states can be seen clearly. 



IV. RATE EQUATIONS 

In the previous paragraphs, we have argued that due 
to the interplay of the spectral distortion and the selec- 
tion rule, there are two time scales present in the prob- 
lem. One describes the tunnel dynamics and is of order 
0{H^). The other is due to higher-order processes and 
is an estimate of how long the system spends in the trap- 
ping states. It is thus evident that the precise implemen- 
tation of the internal dynamics is irrelevant as long as the 
key ingredients to the separation of time scales, namely 
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FIG. 5: Angular-momentum resolved spectra like in Fig. |4l 
In this case the Jahn- Teller molecule is placed slightly below 
the chemical potential of the source electrode. The chosen 
parameters are eV^d ~ 2fitj, eVg = 0.5fkj. Contrary to the 
situation shown in Fig. lU there are no trapping states. 



FIG. 6: Angular momentum resolved spectra for tunneling to 
the drain electrode, eV^d ~ 2huj, eVg — — O.Sfto;. The trapping 
state 1/3) = |1; ^, 0) is marked in red. The three nearest energy 
manifolds of the uncharged system are shaded in gray. The 
dashed arrows show transitions that are forbidden, because 
the target manifold's energy is too large. The dotted arrow 
shows a transition that is forbidden due to the selection rule. 



selection rule and spectral distortion, are included. 

In this paragraph, we present such a numerical im- 
plementation that accounts for the selection rule by em- 
ploying rate equations in the pseudo angular momentum 
basis. By doing so, we neglect all coherences between any 
two states in the density matrix, and therefore also those 
between degenerate states. This will, as we have pointed 
out before, not affect the appearance and structure of 
the current-blockade, because the effect is independent 
of the specific form of the eigenbasis of the stationary 
density matrix. The only error that is introduced by the 
restriction to rate equations is a quantitative one. The 
rate equations, on the contrary, allow for a very intuitive 
and generic treatment of transport problems for systems 
with spectral distortions and selection rules. 



A. Formalism 

The occupation probability P^' of the molecular eigen- 
states |n; q) (n = 0, 1 denotes the charge of the molecule, 
and g is a multiindex for the internal degrees of freedom, 
r is the broadening of the molecular levels due to the 
coupling to the leads.) is governed by the rate equation^ 



5:[ 

n' ,g' 



1 q^q' 



7rol- 



(7) 



The tunnel rates W^^^ are a product of the Golden- 
rule expressions Eq. ^ and the Fermi functions of the 
leads.S'20 The last term, 7,.ei = -(l/r)[P,"-P°^ E,' Pq'l 
describes vibrational relaxation towards the thermal 
equilibrium distribution P^'^ by higher order processes 
not included in the Hamiltonian. By doing so, we ef- 
fectively implement the aforementioned processes, which 
release the system from the trapping states. Our results, 
however, are independent of the details of this implemen- 
tation. The maximal tunneling rate, Tjj~^, defines the 
proper time scale on which we measure all times. 
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The stationary solution of eq. ([7]) permits us to 
compute the stationary current through, say, the left 
junction,— 

= E - , (8) 

qq' 

as a function of the source-drain voltage Vsd and the gate 
voltage Vg. The results for the differential conductance 
dl/dVsd shown in Fig. [3] for (a) fast and (b) slow vi- 
brational relaxation differ strongly from those obtained 
for single-mode models*§^ When the molecular orbitals are 
close to the Fermi energy of the drain electrode (negative 
Vg in Fig. [3]), we find, in accordance with the discussion 
in section [Till an extended region of current suppression 
outside the Coulomb-blockaded region with large nega- 
tive differential conductance (NDC). No such suppres- 
sion occurs when the molecular orbitals are close to the 
chemical potential of the source electrode (positive Vg in 
Fig. [3]), thus rendering dl/dV asymmetric in the gate 
voltage. The positions of the dl / dV-\mes that mark the 
area of the suppressed current in the diagrams are de- 
fined by those values of {Vsd,Vg) for which the trapping 
state is connected to an energy manifold that participates 
in the transport process. If an increase of the bias volt- 
age generates a finite escape rate for population from the 
trapping state, that state is no longer absorbing, and a 
stationary current can flow, hence a line of positive dif- 
ferential conductance is observed. In the opposite case, 
that is trapping behaviour, a line of negative differen- 
tial conductance appears as soon as the principal current 
carrying states — most often it is sufficient to consider the 
vibronic ground state of the system — have tunnel rates 
into the trapping state. Because trapping states are ex- 
cited vibronic states and intermediate excited states are 
necessary to access these, there is a certain energy thresh- 
old below which a finite current may flow but the trap- 
ping state is still out of reach. This threshold will thus 
be visible as a negative differential conductance peaks in 
the dl/dV diagram. Since the trapping states with larger 
pseudo angular momentum have larger energy thresholds 
before they can be populated in steady state, the current- 
blockade mechanism will yield a step-like structure of 
negative differential conductance peaks and plateaus of 
suppressed current. These are clearly visible in our nu- 
merical results in Fig. [3l 

B. Separation of time scales 

We further support the irrelevance of the details of 
the fast, order 0{H^) dynamics by reproducing the po- 
sitions and line shapes of the peaks in the numerically 
obtained differential conductance, see Fig. [31 within a 
simple model. This model only retains the most basic 
features of the dynamics, namely the presence of trap- 
ping states and the two relevant time scales. Consider a 
two-level system with a trapping state \(3) and a state \a) 
from which the electron can tunnel to the drain electrode. 



say the absolute ground state of the charged molecule. 
The current is given by the inverse of the time Ttotai the 
electron spends on the molecule. In our model, this time 
is the sum of the trapping time r and the time the elec- 
tron needs to leave the second, non-trapping state. The 
rate to leave this latter state is the product of the bare 
tunneling rate I/tq and the probability 1 — / to find an 
empty state of appropriate energy in the drain electrode, 
where / is this electrode's Fermi function. The typi- 
cal time for tunneling from \a) out to the leads is thus 
To/(l — /) and the stationary current can be estimated 
as 



Ttotai '^+IT?7) (1-/)t + To' 

This yields the relation dl/dV oc I{tI - l)f/T. Its con- 
sequences for the line-shape are twofold. As a purely 
quantitative effect, the height of the peak in dl/dV de- 
creases with increasing r. As a qualitative effect, how- 
ever, the position of the peak is shifted. In a model 
with just a single time scale, be it a quantum dot or 
a single-mode molecule, the stationary current through 
the drain electrode is only determined by the Fermi func- 
tion 1 — /(ei;q — £o;q' ~ l^a), that is the probability to find 
an empty state in the leads with just the energy neces- 
sary for the vibrational transition \Q\q') ^ ll;?)- The 
maximum of the d//c?t^-peak is thus located at the bias 
voltage for which /ict — £i;q — £o-q' ■ Within our two-state 
model, this behavior is recovered for t = which corre- 
sponds to the absence of any trapping effect due to fast 
relaxation. When trapping is included, and hence r is 
finite, the maximum of the differential conductance peak 
is shifted. 

This simple model compares remarkably well to the nu- 
merical rate-equation results of the previous section. We 
focus on a single peak in the dl/dV diagram and study 
its line shape for different values of the phenomenologi- 
cal relaxation rate r of Eq. ([7|) . The effect of increasing 
T — red to blue curves in Fig. [7| (a) — is a reduction of 
the peak's height and a shift of its position. We can fit 
the results of our simple model, Eq. Q, to the data in 
Fig. [7] (a) , obtaining Fig. [7] (b) . Comparing the fit values 
for r with the parameters of the numerical simulation, 
see Fig. |8l we can see how well this simplified model de- 
scribes the interplay of the two time scales. The only 
deviations occur when both time scales are of compara- 
ble magnitude, r — 0{to). In contrast, for a sufficiently 
large separation of time scales, the agreement is nearly 
perfect. This backs our general considerations of section 
mil where we have argued that the observed phenomena 
are of a generic nature. 

The predominant features of electronic transport 
through a single Jahn- Teller molecule are determined by 
three properties of the system: the spectral distortion 
and the selection rule fix the position of the peaks in the 
differential conductivity, their heights are determined by 
the Franck-Condon matrix elements. Finally, higher or- 
der processes introduce a second, larger time scale, which 
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FIG. 7: Positive differential conduction peak at the Coulomb- 
blockade, Vg = —O.bhu, for different values of the relaxation 
time T (a) and fit with the simplified model, dots in (b). The 
read curves, that is the large, centered peaks have small, the 
blue ones, damped and shifted to lower V^d very large r. 



eventually determines the exact line shape. 



V. EXPERIMENTAL ISSUES 

We finally discuss several experimental issues. Placing 
a Jahn- Teller molecule between two biased electrodes ex- 
poses it to a potentially large electric field. For a semicon- 
ductor quantum dot, weak coupling and good screening 
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FIG. 8: Fit r- values of the simplified model (dots) and its 
application to Fig. [7] vs. the numerical values used in the 
actual simulation. The curve is a linear fit to these values, 
Tfit cx Tnumcricai, showiug the Validity of the simplified ap- 
proach in the regime of strong separation of the time scales. 



implies that the potential drop is concentrated in the con- 
tact regions between the electrodes and the dot.— For a 
single molecule, however, this is less obvious. Thus, the 
electric field may lift the electronic degeneracy by the 
Stark effect and mix the vibronic levels by breaking the 
rotational symmetry of the adiabatic potential. 

To estimate the mixing of vibronic levels, we note that 
the relevant vibrational modes 9 and s do not possess a 
dipole moment and are therefore unaffected by the elec- 
tric field unless anharmonic mode mixing becomes rele- 
vant when other, "dipole-active" modes are strongly dis- 
placed. One readily estimates, however, that for relevant 
electric fields (corresponding to voltage drops of the or- 
der of huj) the typical displacements of dipole-active vi- 
brational modes are small compared with their oscillator 
length ioBc, so that no anharmonic effects are expected. 
Assuming that the dipole moment is of order ex, where x 
is the displacement of the dipole-active mode, we obtain a 
displacement x ~ eE/Muj"^ ~ ^Lc/'^ ^ ^osc- The length 
5 is of the order of the size of the molecular junction. 
We also note that the electric field does not directly cou- 
ple the two electronic levels as both are d orbitals. The 
Stark shift is therefore quadratic in the electric field, and 
we estimate that it is several orders of magnitude smaller 
than the vibrational frequency. 

The coupling of the Jahn- Teller molecule to fermionic 
reservoirs, the electronic leads, may also introduce exter- 
nal perturbations that are able to break the symmetries 
of the system on which our approach relies so heavily. If 
the perturbation is too strong, for example, the degen- 
eracy of the vibrational modes 6 and e is lifted and the 
energies are shifted by a large amount. The Jahn- Teller 
effect is then no longer oi E®e- but rather of i?(8)(5i -1-62)- 
type. A selection rule like the one used in this paper is 
not available in such a model, and the current-blockade 
would not be visible. If, however, the perturbation is only 
small compared with the tunnel coupling, the mixing of 
the states induced by the perturbation is able to intro- 
duce a second time scale t tq into the problem. As we 
have discussed in the previous paragraphs, the presence 
of two widely separated time scales only shifts the posi- 
tions and alters the shape of the differential conductance 
peaks, no matter what the detailed physics behind the 
phenomenology of the conduction peaks are. 



VI. CONCLUSIONS 

We have shown that symmetry-induced degeneracies of 
electronic and vibrational modes can lead to new trans- 
port phenomena in single-molecule junctions. When the 
molecular junction consists of, say, an octahedral com- 
plex with an E^e Jahn- Teller effect, a combination of 
Berry-phase effects and a spectral trapping mechanism 
results in a non-trivial current-blockade, even for param- 
eters where the Coulomb blockade is lifted. It is evident 
from the physics of this blockade mechanism developed 
in this paper that it is associated with strongly enhanced 
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shot noise of the current.^ We have also shown that the 
mechanism is quite generic and that the introduction of a 
second time scale to the problem, for example by higher 
order relaxation processes, co-tunneling, or external per- 
turbations, has only quantitative effects on our results. 

The good agreement of the numerical results of our ef- 
fective rate model with the general consideration of sec- 
tion mil shows that the current-blockade that has been 
found in this context is a generic effect of systems with se- 
lection rules and spectral distortions. Nowack et alJ^ find 
a related current-blockade effect in systems where the 
strong spectral distortion alone damps certain Franck- 
Condon matrix elements. This is, however, no selection 
rule that derives from symmetries like in this paper; it 
rather stems from the almost vanishing overlap of dis- 
placed oscillator wave functions. 



Our discussion of Jahn- Teller molecules raises several 
interesting questions. The most important question con- 
cerns a complete understanding of the role of coherences 
between degenerate states and their implications for the 
transport properties of the system. What are the exper- 
imentally measurable quantities that are influenced by 
non-diagonal elements of the density matrix? And when 
is a reduction to rate equations allowed? 
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